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Abstract 

We discuss a 3D model describing the time evolution of nematic liquid crystals in the framework 
of Landau-de Gennes theory, where the natural physical constraints are enforced by a singular free 
energy bulk potential proposed by J.M. Ball and A. Majumdar. The thermal effects are present 
through the component of the free energy that accounts for intermolecular interactions. The model is 
consistent with the general principle of thermodynamics and mathematically tractable. We identify 
the a priori estimates for the associated system of evolutionary partial differential equations and 
construct global-in-time weak solutions for arbitrary physically relevant initial data. 
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1 Introduction 

The main aim of this paper is to derive and analyze a thermodynamically consistent system of evolu- 
tionary equations describing the dynamics of nematic liquid crystal flows in 3D. We use the abstract 
thermodynamic framework proposed by Frcmond |12| in conjunction with Beris-Edwards formulation [3] 
of isothermal liquid crystal hydrodynamics (cf., e.g., |7]). The state of the complex fluid at the time t 
and spatial position x is described by means of a Q-tensor field Q = for the nematic director 

orientation, the velocity field u = u(t,x), and the absolute temperature § = i?(t,a;). 

The main characteristic of nematic liquid crystals is the locally preferred orientation of the nematic 
molecule directors. This can be described by the Q-tensors, that are suitably normalized second order 
moments of the probability distribution function of the molecules. More precisely if is a probability 
measure on the unit sphere S^, representing the orientation of the molecules at a point x in space, then 
a Q-tensor (Q)(a;) is a symmetric and traceless 3x3 matrix defined as 

and it is supposed to be a a crude measure of how the probability measure fix deviates from the isotropic 
measure p, where dft = ^dA, see [6]. In the Onsager model (cf. [6], [T9|) the probability measure is 
assumed to be continuous with density p = p(p). Then 

Q{x)^ l^^(p®p~h^ p{p)dp. (1.1) 
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The fact that Hx is a probabihty measure imposes a constraint on the eigenvalues of Q, namely that 
they are bound between the values — ^ and |, see p]- Thus not any traceless 3x3 matrix is a physical 
Q-tensor but only those whose eigenvalues are in (— i, |). The hydrodynamic models available in the 
literature do not have, to our knowledge, a natural way of preserving this physical eigenvalue constraint 
on the traceless and symmetric matrices. One possibility is to use a singular potential, such as the one 
proposed by Ball and Majumdar [5], that enforces the physical constraints satisfied by the Q-tensors, and 
this is the solution we adopt in our model. 

The hydrodynamic theory for the Q-tensorial isothermal model in case of regular bulk potential 

MQ) = ^tr(Q2) _ ^tr(Q-^) + ^tr'{Q') 

has been recently studied in [22] and |23j . 

Let us notice that in the literature there are few papers dealing with non-isothermal models for liquid 
crystal dynamics. Two attempts, in the case when the evolution of the director is described by the 
vectorial director field d (standing for the preferred orientation of the molecules at any point), were 
made in [TI] and [5]. In particular, in [5] the stretching and rotation effects of the director field induced 
by the straining of the fluid were considered and the existence of global in time weak solutions was 
obtained for the corresponding initial boundary value problem. In the present contribution we follow 
the thermodynamic approach exploited in [8] in order to deal with the tensorial model obtained using 
a non- isothermal version of the singular bulk potential /(Q) proposed in [2]. In the spirit of [6], we 
include the temperature dependence in the potential assuming that the coupling term in the free energy 
functional is given by (cf. also [25] and [20] ) 

V'b(i?,Q) = ./(Q)-C/WG(Q). (1.2) 

Hence, ipB is the sum of a singular part /, independent of temperature with a smooth perturbation 
depending both on i) and on Q. We assume U to be a convex and decreasing function of i!) having controlled 
growth at infinity (cf. (|1.6H1.7|) below). Actually, in a neighbourhood of a characteristic temperature i?* 
of the crystal the function U can display a linear growth and can change sign at "d*, like U {'&) = a{d* —"d). 
According to [6], the function G can be, e.g., given by G(Q) = Q^Qij = tr(Q^). Let us notice that in 
the Ball and Majumdar [2] approach the temperature dependence is of a different type: actually, they 
assume 

M^,Q)^^fm + Gm- (1-3) 

However in [2] only the stationary case is considered; hence, one can freely divide by ■& the expression 
tpB in (|1.3p and obtain an expression which corresponds to (|1.2p at least for values of not too distant 
from the critical temperature It is worth noting that, in the evolutionary setting, dealing with a 
free energy of the form (|1.3p would be mathematically more complicated since the coupling occurs in the 
singular part of the potential; for this reason we expect that weaker analytic results could be proved in 
that case. We will devote a forthcoming paper to the analysis of the evolutive model with the free energy 

Comparing the present analysis with the previous paper [8], a major difficulty is provided here by 
the presence of the singular potential /, which has to be properly handled by means of convex analysis 
tools. In addition to that, we consider here a more complicated version of the heat equation involving an 
explicit dependence of the thermal relaxation coefficient with respect to the Q-tensor (cf. (|4.7p below). 
Actually, this choice, which is more realistic from the physical point of view (in particular, it gives rise 
to an entropy s depending also on Q and not only on d, cf. ()1.40p below), creates a number of additional 
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mathematical difficulties. The key point, which requires some care to be accomplished, is related to the 
proof of regularity and strict positivity of d at the approximate level, two properties which are crucial 
for the purposes of proving the validity of the entropy inequality and of the total energy balance in the 
frame of weak solutions. 



1.1 Landau-de Gennes free energy with the Ball-Majumdar bulk potential 

Denote by Rgy^^Q the linear space of symmetric traceless 3x3 real- valued matrices. The Landau-de 
Gennes free energy takes the form 

^(Q, V,Q, I?) = i|V,Qp + i^Bi-d, Q) ~^log{d), (1.4) 

where Q{x) e for aU X in the smooth domain SI C R . 

Ball and Majumdar [2] introduced the bulk component of the internal energy functional by means of 
a singular functional t/'b = ipsi'^, Q) that, for any fixed temperature ??, blows up when at least one of the 
eigenvalues of Q approaches the limiting value —1/3. In particular, the boundedness of the free energy 
enforces the boundedness of Q in L°° . Specifically, we set 



where 



V'b(^,Q) = /(Q) - U{d)Gm for Q e Rl^l^, (1.5) 

infpe^^ p(p) log(p(p)) dp if KM e (-1/3, 2/3), i - 1, 2, 3, 
oo otherwise. 



p G L\S^), J^^ p{p) dp = 1; Q = ^ (^p ® p - p(p) dp 



The function / is the singular component of the bulk potential. In here singular refers to the fact 
that the domain is not the whole space (while inside the domain the function is in fact smooth). The 
function / enjoys the following properties that can easily be deduced from [21 Section 3, Prop. 1]: 

• / : i?s^„i ^ [^^7 °°] is convex and lower semi-continuous, with K > 0. 

• The domain of /, 



V[f] = {Q G i?f;„to I /(Q) < ^} = {Q G Rt^l, I \M G (-1/3, 2/3)}, 

i3x3 
'sym,0 ■ 



is an open convex subset of R'^^^ 



• / is smooth in 

The potential G characterizes the action of intermolecular forces. In contrast with Ball and Majumdar 
[2] , we suppose the temperature changes act on this component of the bulk potential. Such a hypothesis 
is quite common in the literature (see for instance [5]); actually, U is typically assumed to change sign at 
a critical temperature. Here, we assume that there exists a positive constant c such that 

U eC°[0,+oo)nC'\0,+oo), C/(0)>0, U' <0, [/ convex in [0,oo), limsup C/"(?9)t93/2 < +oo, (--^ g^ 

1?— >-oo 
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\U\d)\<c\§\-^/'^ for all i9 e (0,oo), (1.7) 

and 

G e C3(i?3x3^)^ G > 0, G(Q) = G{RQR') for all R e SO{3). (1.8) 

Let us note that the choice G(Q) = tr((Q^) corresponds to that of 0. Moreover, it is not restricive to 
assume that G is uniformly bounded together with his first and second derivatives. 

1.2 Thermodynamics 

In accordance with the general thermodynamical framework of |12] , we introduce the set of state variables 

^ = (Q,V,Q,i?), 

together with the dissipative variables 



where 

^(u) = ^(V.u + V^u) 

is the symmetric velocity gradient, and 



= a^Q + u • V,Q - §(V,u, Q) (1.9) 
is an analogue of material derivative characterizing the time evolution of the tensor Q, with 

§(V,u, Q) = (e£(u) + c^(u)) + il^ + - u{n)) - 2^ (^Q + ^I^ (Q : V,u) , (1.10) 

w(u) = ^(V.u-V^u), 

where ^ is a fixed scalar parameter, measuring the ratio between the rotation and the aligning effect that 
a shear flow would have over the directors, see Beris and Edwards [3]. 
The evolution is ruled by the pseudopotential of dissipation $ 



<i>(£;^£;)=/.(^)|£(u)|2+/o(div,u) + 4?|V,^^P ' ^ 



' ' 2r(i9) 



Dt 



(1.11) 



with the shear viscosity coefficient the heat conductivity coefficient k, and the collective rotational 
viscosity coefficient T. The incompressibility of the fluid is formally enforced by the functional Iq - the 
indicator function of the point {0}, 



loiz) 



if z = 0, 
+CXD otherwise. 
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Similarly to U, the transport coefficients ^J.{i}), K{d), and T{'d) cliange with temperature. For the sake 
of simplicity, we suppose that 



jj., K, T e C^[0, oo), < 



' < ^ < l-iW < 
< K < < K 

. < r < r(i9) < r 



\ for all i9 > 0. 



(1.12) 



1.3 Time evolution 



We assume that the fluid has a constant density g, say g ~ 1. Then, in accordance with the general prin- 
ciples developed in Fremond [12l Chapters 2,3], the time evolution of the system is uniquely determined 
by the choice of the potentials T and 3>. 



1.3.1 Momentum equation 

Newton's second law is expressed by means of a modified Navier-Stokes system: 

dtu + div2;(u u) = divxO- + g, 
where g is a driving force, and a denotes the stress tensor. 



^ I nd 

(7 — a + a 



The dissipative component of the stress reads 

d 5$ ^(^9) 



de{u) 



(V^u + Vlu) - pi, 



where p is the pressure. We have 



-p e 5/o(div2;u) 
yielding the standard incompressibility constraint 

div^jU ~ 0. 

The specific form of cr"'^ will be derived below. 

1.3.2 Entropy production 

The heat flux q can be decomposed as 

q = q' + q"^ 

where the dissipative component obeys the standard Fourier law 



(1.13) 



(1.14) 



with the associated entropy flux qe = /■&. 
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The energy density tensor is taken to be 



f^ pQ Dt 

^ Dt 



where 



piid 



-tr 



c 



U{d)C 



(1.15) 



for any 



e R. 



3x3 
sym ' 



denotes the projection onto the space of traceless tensors. In other words, B"*^ can be seen as the 
subdifferential of J- with respect to Q in that space. 

We state the entropy equation in the form (cf. [5J [5] for a complete derivation of this equation) 



dts + div,(5u) - div.qe = ^ : + ®^ • ^ + ^l^-^P 



(1.16) 



Dt 



V,i?p >0, 



with the entropy 



l+l0g(79) + [/'(79)G(Q). 



(1.17) 



Note that, in accordance with hypotheses ()1.6p . (jl.Sp . the entropy s is an increasing function of the 
temperature. 



1.3.3 Q-tensor evolution 

The internal energy balance reads 



dte + div2,(eu) + divj;q = cr : Vi:U 



Bt ^-^--^^'Dt 



where 



e = P+ = -|V,Q|2 + /(Q) - - m'i§) 



and Y is the energy flux tensor 



Y = Y"'^ 



(1.18) 
(1.19) 
(1.20) 



For simplicity, we assume here that Y has no dissipative component (and, correspondingly, that $ is 
independent of Va;Qt). 

The principle of virtual powers (see Fremond jl2l Chapter 2]) yields the time evolution of Q, namely, 

div^Y = B; (1.21) 



in other words, 

dtQ + div,(Qu) - S(V,u, Q) = r(i9) [AQ-C 



(1.22) 



where S is defined in (|1.10p . It is easy to check that the space Rly^ q is invariant for solutions of (|1.22p . 
specifically Q(t, •) G i?sy1n.o fo^' t > as soon as Q(0, •) € R\ 



3x3 
sy m , ■ 
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1.3.4 Total energy balance 

Taking the scalar product of the momentum equation (|1.13p with u and adding the resulting expression 
to (jl.lSp . we obtain the total energy balance in the form 



dt -|uP + e +div. 



up + e u + div^q 



(1.23) 



diva,(cru) + div^j T{d)VxQ ■ 



gu. 



Dt 



(1.24) 



(1.25) 



It remains to determine q'^'^ and cr"'^. 

To this end, wc first multiply the entropy balance (|1.16p by -d. This gives 



-dst - su • + div^(t9su) + div.T;q'' = a'^ : e(u) + I 

Next, using (dUT]), (ITigi) . p:T5)) and ([TT^ . we get 

et=Y::atV,Q + B"'i + 
Moreover, 

div,(eu) = div,((J- + st?)u) = u • (J-tfV.T? + J-Q : V,(Q + J^v^q :: V,V,Q) + div,(si9u) (1.26) 

= -su • V^d + u • (B"'^ : V^Q) + u • (Y :: V^^^V^^.Q) + div^,(s'i9u). 
Replacing (|L25Hr26l) in ([TTTS)) and using ((L9)) . wc then have 

Y :: StV^Q + B"'^ : (9tQ + 7?st - su • V^i? + u • (B"'^ : V^Q) + u-(Y :: V^V^^Q) +div:r(si?u) + div^q (1.27) 



= a'^ : V:,u + cr"'^ : V^u 



(7"" : V:rU + cr"'' : V^u 



Dt Dt 

nd 



+ B"'^ : + u • V,Q - §(V,u, Q)) + Y :: (9t V^Q + V,(u • V,Q) - V,§(V,u, 
Simplifying some terms and using symmetry of a'^ , we have more precisely 

■dst - us • Vx-d + divj,(st?u) + div^q'^ + div^q"'' 



(1.28) 



= cr'^ : e(u) + ct"'^ : V^^^u + 1"^ : — ^ - B""* : ^(V^rU, Q) + Y :: (V^u ■ V^^Q) - Y :: V:r§(V^u, ( 
Then, subtracting (|1.24p from (|1.28p and using (|1.15p . we arrive at 

div,q"d 



(1.29) 



a"^ : V,.u - ( C 



- U{d)C 



S(V,u, 



V:,S(V:,U, Q) + (V:,Q V^Q) : V:,U 



= cr"'^ : V,.u - div^ {V^-Q : §(V^u, 



-C 



U{d)C 



S(V^u, Q) + (V^Q V^Q) : V^u. 



Consequently, we deduce that 
and 



-V,Q:§(V,u,( 



HQ + 2e [H : Q] ( Q + -I ) - e 



where we have denoted 



and we have used the identity 



U{i))C 



(1.30) 

(V,Q0V,Q), (1.31) 

(1.32) 
(1.33) 



3 



: V ;^u 



-H: §(V^u, 
V,u + 2^ (H : Q) (Q : V,u) - ^ 
that holds for any symmetric matrix H. 
1.3.5 Evolutionary system 

The computations given in the previous section permit to write the resulting evolutionary system in a 
concise form: 



Incompressibility: 
diva;U = 0; 

Momentum equation: 
9tU + divj; (u u) = div^cr + g; 

Order parameter evolution: 
+ cliv,(Qu) - S(V,u, Q) = T{-d)m- 



dt 1 ^\n[ 



Total energy balance: 
'1 



div.T 



u ■ 



div^(cru) +div^(r(z?)V:r 



e u + divj^q 



(1.34) 
(1.35) 
(1.36) 

(1.37) 



with the stress tensor 



a = n{^) (V^u + V^u) - pi 



(1.38) 
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+2e(H:Q) (Q+-I) 



where 



L 



'0 V:, 



and the internal energy flux 
where 



-/t(i?)V,i9-V,Q:§(V,u,( 



(1.39) 



S(V:,U, 



1 



V:,U) 



We recall that 



Ui-d) - i9C/'(i?)) G(Q) +1?, s = 1 + log(i9) + U'{^)G{Q), (1.40) 



where we have anticipated that the relation C[(J] = Q is preserved in the course of evolution. 
The system p.34ffTT37|) may be supplemented by the entropy inequality 



> 



dts + div2;(su) — diVa: 

1 ffii^) 



k(i9) 



(1.41) 



|V,u + V^u| +r(i?) |H|- 



In order to avoid problems related to the presence of a kinematic boundary, we suppose that the fluid 
motion is spatially periodic. This can be conveniently formulated by taking the spatial domain f2 C M"^ 
as a flat torus 



0= ([-^,7r]|{_,,,})'. (1.42) 



Note also that the pressure p appears explicitly in the energy balance (|1.37p . in particular, p must be 
determined from the Navier-Stokes system (|1.35|) by means of the Helmholtz projection. Such a step may 
involve insurmountable difficulties in the case of general boundary conditions. 
The original state of the system is given by the initial conditions 



u(0, •) = uo, (Q(0, •) = Qo, ^(0, •) = ^?o. (1.43) 



Our goal in this paper is to study the initial- value problem (jl.34ffT37|) . supplemented with the bound- 
ary conditions (|1.42[) and the initial conditions (|1.43p . in the framework of weak solutions. In the sequel 
we will show that, for any choice of finite energy initial data (cf. Section [2] below), the problem possesses 
a global-in-time weak solution, which additionally satisfies the entropy inequality ()1.4ip in the sense of 
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distributions. To this end, we first derive formal a priori bounds in order to facilitate the reading of the 
main rather technical part of the proof, see Section [3] The global-in-time weak solutions are constructed 
as a limit of solutions of a family of approximate problems introduced in Section |4l The most delicate 
part of the proof is showing strict positivity of the (absolute) temperature by means of a weak variant 
of the parabolic comparison theorem. The proof of convergence of approximate solutions is completed in 
Section O 

2 Weak solutions, main results 



Weak solutions to the problem 34141. 37|) . ()1.42|) . ()1.43p belong to the regularity classes indicated by 
the a priori bounds discussed in Section [5] below. In particular, we have 



u e L°°(0,r;L2(17;M3)), < 



QeL°°((0,r)xr!;i?3;<^%), 
V,QeL°°(0,r;L2(r!;i?27)), 



}, I z?eC^([0,r];i/-2(i7)) 

_ log{d)eL°-{0,T;L\n)) J 



y (2.1) 



and 



V,ueL2((o,T) xl];i?3x3), QeL'{0,T;W''\n;Rl^l,)), V,^? G L«((0, T) x 17; i?^) (2.2) 
for any q < 5/3. 

2.1 Weak solutions 

The weak solutions are defined in the standard way. Given the anticipated regularity of the velocity 
field, the incompressibility constraint (ll.34p makes sense a.e. in the set (0,r) x fl, while the momentum 
balance (I1.35P is replaced by a family of integral identities 







u • dtip -I- (u (g) u) : 



dx dt 



a : Vi,(/7 — g ■ </?) dx dt — / Uq • (p{0, •) dx (2.3) 
"'f2 Jn 



satisfied for any test function ip e C^{[0,T) x Q;R^) and where a is defined in (|1.38p . Here and in 
what follows, we always tacitly assume that all quantities appearing under the integrals are (at least) 
summable in (0,T) x fi. 

Similarly, the evolutionary equation ()1.36p for the Q-tensor is replaced by 



T 

Jn 



dtf+[uQ] : V:,^ + §(V^u, 



dx dt= - / r{i})M -.ipdxdt- / Qo : vi^, •) dx 
Jo Jn Jn 

(2.4) 

for any ip € C^{[0,T) x f7;i?3x3)^ ^y^^^^ S{VxU,Q) is defined in (fLlO)) and H in ![L32^ . 
The total energy balance (|1.37p is satisfied in the sense of integral identity 



T 

Jn 



dx dt 



(2.5) 
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for any ip S C^([0, T) x fl), where q is defined in (|1.39[) and e in (|1.40p and we have set 
eo - ^iV.Qol' + /(Qo) - - ^oU'i^o))G{Qo) + ^o- 

2.2 Main result 

We are ready to state the main result of this paper. 



Theorem 2.1. Let the initial data Uq, Qqj '^'^'^ '^o be given such that 

uo G L^(r2;i?^), div^uo = 0, < , z?o S L°°(f^), essinf t9o = j? > 0. (2.6) 

Suppose that the functions U , G satisfy the hypotheses i [J.61jJ.7[ ), IJ.ig)) . the transport coefficients /z, 
K, and r comply with \1.12]) . and that 

geL"-{0,T;L'{n;R^)). 

Then the problem Ijl.SjWl.ST^ , {1-4^ , il43\j admits a weak solution u, Q, -d in (0, T) x in the 
sense specified in Section \2.1[ In addition, there exist positive constants c and A such that 



I ■d{t,-)>cexp{-Xt)dfor aUt>0, (2.7) 

and the entropy inequality |-?.^j[ ) holds in the sense of distributions. 



The rest of the paper is devoted to the proof of Theorem 12.11 For the sake of simplicity, we set g = 
as the proof in the more general case requires only straightforward modifications. 

3 A priori bounds 

A priori bounds are natural (formal) constraints imposed on hypothetical smooth solutions by the equa- 
tions and the initial data. 

3.1 Energy bounds 

Uniform boundedness in time of the total energy is straightforward consequence of (jl.37|) . In accordance 
with ([2T|) . we get 

uG L°°(0,r;L2(f7;i?3)), (3.1) 
e L°°{0,T-,L^{n)), in particular, Q e L°^{{0,T)xn;Rl^^ o)^ ya^QeL°°{0,T;L^{n;R^^)), (3.2) 
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and 

I? e L°°{0,T;L\n)), 

where we have anticipated the fact that the absolute temperature is a positive quantity. 

3.2 Entropy bounds 

Integrating the entropy inequahty (|1.4ip and using p.3p . (|1.12p f61 we infer that 

\og{-d) e L°°(0, T; L\n)) D L^{0, T; W^^^{n)). 

3.3 Bounds based on energy dissipation 

Muhiplying the entropy inequahty (jl.4ip by ?? we deduce the thermal energy balance in the form 



(3.3) 



(3.4) 



dti^ + u- - d 



i9 > -i9 



dt(U'i^)GiQ)) + u . VJ [/'(i?)G(Q) 



M(^) 



^ |V:rU + V>| +r(i?)H:H, 
where, furthermore, 

^dt(u'{'d)G{Q)) = dt [[m'i^) - U{d) + (7(0)) G(Q)] + (U{d) - C/(0))atG(Q), 



whence 



dt (z9C/'(i9)-[/(i?) + [/(0))g(Q)+i9 + u • 
liv,(K(i9)V,i?) > (c/(0) - f/(i?)) [atG(Q) + u • V, 

= (uio)-uid))c 



(t9[/'(i9)-[/(i9) + t/(0))G((Q)+z9 (3.5) 
m(^?) 



V^u + V> 



S(V,u,Q) + r(i9)H 



V^u + Vlu 



r(?9)H : H 
r(z9)H : H. 



Since we aheady know that Q is uniformly bounded, we deduce from (j3.5p that 



> 



dt 
m(^) 



(?9C/'(i9) - ?/(i?) + C/(0))g(Q) + ?9 



da; 



(3.6) 



c 



U{0) - U{^) 



where G > is an explicitly computable constant depending only on ||Q||l°°, G and ^. 

By virtue of hypotheses (ll.6H1.7p . the function U'{'d) is bounded, therefore we may combine 
with the energy estimate p.3p . use assumptions (|1.12p and the convexity of / to conclude that 



V,ueL2((0,r) X 17;i?3), QeL^{Q,T;W^'\n;Rl:'l^)), and/: 



dfm 



eL\Q,T-L\n-Rl^l^)). 



Finally, thanks to (|3.H43.2p and p.7p , a comparison of terms in (jl.36p gives also 

dt^^L\0,T-L%a-Rl^l^)). 



(3.7) 
(3.8) 
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3.4 Bounds on the temperature gradient 

The estimates on the temperature gradient are obtained by multiplying (|3.5p on —(1 + 19)"", a > 0. Note 
that 

Observe also that the right hand side of p.Sp is uniformly bounded in thanks to the previous estimates. 
Moreover, as we integrate by parts the terms depending on U, we obtain, on the right hand side, the 
terms 



which we need to control (as well as similar quantities depending on the transport part of the material 
derivative, which can be treated in the same way). Integrating by parts in time, we get 



dt{G{q)Y^{d)) 

where 



(1 + ^) 



i3 



G(Q)t, (3.9) 



YM) = / :-U"(s) ds 

^ ' Ji (1 + s)" ^ ' 



and, thanks to (|1.6H1.7p . the function in square brackets in (13. 9p goes like -d^^"^ " for large d. Hence, due 
to dill), it lies (at least) in L°°{0,T; L'^{n)). 

Now, using that G is bounded with its first derivatives and recalling <\3.8\i . all terms in (|3.9p can be 
controlled. Hence, after a straightforward manipulation we can conclude that 

Va;(l +1?)^ e i^((0,T) X n;R^) for any a > 0. (3.10) 



The a-priori bounds obtained here are enough to make the weak formulation (j2.3H2.5|) meaningful. 
In particular, the pressure p (cf. (|1.38|) ) can be "computed" directly from (|2.3p (cf. 9 and 8^ for more 
details) and it is possible to obtain 

L5/^((o,r) X n). 

It can be shown that the a priori bounds obtained in this section are strong enough in order to 
establish the weak sequential stability of the family of solutions to our problem. However, we do not 
pursue this path and pass directly to the construction of a family of approximate solutions. 

4 Approximate problems 

The weak solution, the existence of which is claimed in Theorem 12.11 will be constructed by means of a 
family of approximate problems. 

4.1 Approximate velocity fields 

The velocity field u is obtained via the standard Faedo-Galerkin method based on the family of finite 
dimensional spaces 

Xn = |v £ C°°(fl; R^) V — a trigonometric polynomial of order N, div^^v — o| . 
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Accordingly, the momentum equation (|2.3p is replaced by a finite system of ordinary differential 
equations 

— / u • V dx (4.1) 
Jn 

= / [u]s ® u : V^w Ax- 5 I |V:rUr"^Va;U : Vj;V dx- ^(i?)(Vj;U + V^u) : V^w dx 
Jn Jn Jn 

+ I (V,(Q V,q) : V,v 



dx 



2C(M™,5 :Q) 



with the initial condition 



for any v £ Xat, where 



rii ^^m,> 



/ u(0, ■) • V dx = / [uo]^ ■ V dx 
Jn Jn 



Va;V dx, 



(4.2) 



df„ 



dGiQ) 



(4.3) 



Here r e (3, 10/3), S and m stand for positive parameters; [uj^ denotes the standard regularization 
with respect to the x- variable by means of a family of convolutions, while {/m}m>o is a family of smooth 



convex functions defined on R^y^ g such that 



fm<f foraUmeN, /^^(Q) </„,(< 

fm ^ /, uniformly on compact subsets of T>[f] 



for all mi < TO2, for all Q € Rf"^ 



dfrr 



c 



df_ 



uniformly on compact subsets of 



'sym,0 ' 



> as TO — > OO, 



/„, ^ OO uniformly in R^^^^ q \ 



ci < 



< Cl\Q\ + Cl for all Q e i?f xs^q, to > 



Notice that the term (5| Vajul*" ^ V^u guarantees additional regularity of u needed in the Q-tensor equation 
(cf. estimate (I4.14p and Section [5751 below) . 



4.2 Q-tensors 

The equation governing the time evolution of the approximate Q-tensors reads 

"a/„(Q 



dtQ + (u • V,)Q - §(V,u, Q) = [Vi-d)], (aQ - 



c 



with 



Ms- 



(4.4) 
(4.5) 
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Here, [r(i?)]e denotes a regularization (via convolutions) of T{'d) with respect to both t and x variables. 
Furthermore, [QoU € C°°(f^; i?fy^ g), such that 

[Qo]* ^ Qo a.a. in n, fHQoh) ~^ /(Qo) in L\n). (4.6) 

An explicit construction of [Qo]<5 can be obtained, for instance, by truncation and mollification (recall 
that the domain of / is an open set). Then, the convergence property in (j4.6p can be verified using the 
dominated convergence theorem. 

Finally, Us : M ^ K is a bounded truncation of U satisfying ([LMTTI) and such that Ugi'd) = Ug{0) 
for 1? < 0. 

4.3 Thermal energy balance 

The approximate temperature is determined via a "heat" equation of the form 

dt^ + u- - div^ (k(i9) V:ri9) (4.7) 



\7^U + Vlu 



2 

+ [r(i?)]eH„,5 : H„,5 + ,5|V,ur, 
where IHIm^^ is defined in (|4.3p . Equation (|4.7p is complemented with the initial condition 

z9(0,-) = [^9o]5, (4.8) 
where, similarly to (j4.2p . [ffols denotes a regularization in the space variables. 



4.4 Existence of approximate solutions and uniform bounds 

Our program for the remaining part of the paper will be to construct approximate solutions to the 
problem (|4.1H4.8p and let successively 

m oo, N oo, e 0, and, finally, J — 

in order to recover in the limit a weak solution to the problem (|1.34| - [T757)) . (|1.42| - nn5)) . the existence of 
which is claimed in Theorem 12. II 

For fixed values of the parameters m, N, e, and 6, we can construct local-in-time solutions to the 
approximate system by means of a Schauder fixed-point argument. This procedure is similar to the 
one sketched in [HI Section 5] (see also |TOl Chapter 6] for further details); hence we leave it to the 
reader. Moreover, the local solutions can be extended to the whole time interval [0, T] as soon as suitable 
uniform estimates analogous to the a priori bounds are established. Actually, to simplify notations, we 
shall directly assume that solutions are defined on the whole (0,r) already in the approximation. 
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4.4.1 Energy bounds 

Bounds on the total energy are obtained in the same way as the a priori bounds. We take v — u{t, •) as 
a test function in (|4.ip . multiply (|4.4p by — H^^^, and add the resulting expression to (|4.7p to obtain 



d 



We deduce that 



ess sup \\u{t, •)llL2(n;_R3) < c, 
te(o,T) 



dx = 0. 



ess sup 

te(o,T) 



and 



ess sup \\-dit,-)\\L^n) < c. 
te(o,T) 



(4.9) 

(4.10) 
(4.11) 

(4.12) 



In addition, since u ranges in the finite dimensional space Xn consisting of smooth functions, we get 

sup ||u(i, •)||c'=(0;_R3) < c(fc,iV) for any fc = 0, 1,. . . (4.13) 

te[o,T] 

4.4.2 Bounds on the Q-tensors 

Since the approximate Q-tensors satisfy equation (|4.4p . where the leading coefBcient [r(t?)]j is smooth 
and the coupling term u ■ V^Q — §(V2;U,(Q) is regular thanks to (|4.13|) . we may bootstrap the maximal 
regularity estimates of L'-type (see e.g. Krylov [T5]) to deduce that 



II^*<QIIli(0,T;Li(O;_r3x_3 J) + 



L9(o,T;W2.9(n;_R^;<^^^)) ^ c{q,N,m,e,S) for any 1 < q < oo. (4.14) 



Now, we can go back to (|4.ip and to use (j4.14p to conclude that 

ess sup ||9tu(t, •)||cfc(0;_R3) < c(fc, A^,m,e,(5). 
te(o,T) 



(4.15) 



4.4.3 Strict positivity of the absolute temperature 

This is one of the most delicate steps in the proof of Theorem 12.11 Our aim is to prove that d is strictly 
positive already at the approximate level and uniformly with respect to all approximation parameters. 
To this aim, we apply the parabolic comparison theorem to equation (j4.7p written in the form 



dG{Q) 



dt^ + U-V^^) -div,(K(l?)V,1? 



(4.16) 



(§(V,u,Q) + [r(i9)]eH™.5 



+ [r(i?)]eH„M : H„,^5 + (5|V,ur. 

Now, a short inspection of (jl.lOp yields that all terms can be expressed as products of scalar quantities 
depending on Q and the symmetric gradient e(u) with the only exception of the commutator 
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Fortunately, thanks to assumption (jl.Sp . we have 



u{u)Q - Qcj(u) 



2 C 



dGiQ) 



w(u) = 0. 



Let us now observe that, by prMTTT]) . there exists c > such that |C/'(i?)?9| < cd^/'^ for all d > 0, and 
the same inequality can be required to hold for the approximations Us- Moreover, we can suppose G 
bounded because the domain of / is bounded and consequently the behavior of G for large |Q| is not 
relevant. On account of these considerations, it is not difficult to arrive at 



(l + z9C/ni9)G(Q)) (dt^ + u . V,z9) - div,(K(i?)V,i?) 



(4.17) 



> -Ai9 + 



W^u + Vlu 



[m] 



where A is a positive constant depending only on the structural properties of the functions Us and G. 

Testing (|4.16p by i?" and using the fact that Ug = for < 0, we obtain, by standard arguments, 
the non-negativity of 'd. 

Next, multiplying equation (I4.16P by n{'d)dt'& and making use of the available bounds (I4.12p . (|4.13p . 
(|4?T4l) . and ([LT]) . we deduce 



ess sup ||i?||vvi,2(o) < c(iV,TO,£,(5), \\dt-&\\L-^({^a,T)xn) < c{N,m,e,5). 



(4.18) 



Thus, introducing 



e = /C(?9), /C' = K, /C(0) = 0, 



and noting that the coefficient (1 + -dl/'g' {^)G{Q)) is uniformly bounded, applying standard elliptic regu- 
larity results to (|4.16p we arrive at 



I|0||l2(o,t;W2.2(o)) < c{N,m,e,S). 
whence, by (j4.18p and interpolation, it is not difficult to obtain 

II^IIl2(o,t;W2-3/2(o)) < c{N,m,e,6). 
Next, we may rewrite (|4.17p in the form 

a(t, x) (dtO + xe + u- v^e) - Ae = 5 > 0, 

where the previous estimates imply that g G ^''((0, T) x il) for all p G [1, oo) and we have set 

_ l + m's'{^)G{Q) 



(4.19) 



(4.20) 



(4.21) 



GL°°((0,T) xl7), a>a>0, 



where A > depends only on A and the structural properties of the function k, and a depends only on 
K. For instance, one may take 
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Then, multiplying (|4.2ip by e'*'* and setting z = e'*'*0, we get 



a(i, x) [zt^'\x■ V:rzj - Az = ge^* > 0, (4.23) 

whence, testing by — (z — 8)~, where = /C(?9), and integrating by parts the terms with a, we obtain 
the desired conclusion 

d{t,-) > cexp{~\t)^ for ah t G [0,T]. (4.24) 

Let us note that in principle estimate (|4.24p holds for smooth a (see e.g. [13]). However, considering the 
class of solutions specified through (|4.18p . (|4.20p . the procedure can easily be extended to any bounded 
measurable a by means of density arguments. Indeed, we can consider a smooth (Holder continuous) 
approximation of a and a smooth approximation of g, and we can assume ah ^ a, gh ^ 9 in 
for every p G (1,(X)). Then, by the standard maximum principle argument (cf. |13| ) we can conclude 
that there exists 'dh € H^{0, T; L'^iQ)) n L'^{0, T; W'^^^{n)) (in fact, 'dh wih be even smoother) satisfying 
(|4.24p . Then, also the limit of 'dh satisfies the same inequality. This concludes the proof because we 
have uniqueness of solutions for (|4.2ip . which can be obtained simply testing the difference of the two 
equations by the time derivative of the difference of two solutions. An analogous proof of the maximum 
principle argument for (|4.21[) with bounded measurable a can be found, e.g., in 14, Prop. 3.6, p. 293]. 

Estimate (I4.24p coincides with (|2.7p claimed in Theorem 12.11 It is remarkable that (|4.24|) is indepen- 
dent of the parameters to, N, e, and 6. Indeed, the choice (|4.22|) of A is independent of all approximations. 



4.4.4 Moser and regularity estimates on the absolute temperature 



We prove here that the absolute temperature is (globally in time) Holder continuous at least at the 
approximate level. Although this fact is essentially a consequence of well-known techniques for parabolic 
equations (see, e.g., [E]), in view of the fact that (|4.16|) depends on in a somehow intricated way we 
give, for the convenience of the reader, at least the highlights of a direct proof based on Moser iterations. 
First of all, we rewrite (I4.16P in the form 



1 



m's'id)G{Q)') [dt^ + u • V^i?) - div, (k(^?) V.T?) = ~m's{t9)£ + y, 



(4.25) 



where the functions i and v collect the various quantities on the right hand side of (j4.16p and satisfy 

¥\\li{o,t-m{^)) + \W\\Li{o.T-Li{ii)) ^ c{q,N,m,e,S) for any I < q < oo. (4.26) 

Then, by the lower bound (|4.24p we are allowed to test ()4.26l) by 1)^^^ for a generic p > I. Integrating 
with respect to space variables, we then obtain 



dt 



P 



(4.27) 



{dtQ + u-V^ 



dx. 



where we have set 



Hs^pi^) := / rPU^'ir) dr 
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and we can notice that, due to (|1.6H1.7p (or, more precisely, the analogue for Us), we have 

Then, owing to the fact that Hs^p{{})G{Q) is nonnegative, and using the regularity of Q and u to estimate 
the right hand side, we see that (j4.27|) assumes in fact the structure 

5t||^||^.(o)+ci||V.^^/2||i.(^)^^3 <C2(l+p)^ (0^1/2) (1 + ^^-1) dx, (4.28) 

where ci , C2 > are absolute constants and cj) collects all the terms depending on u and Q (also through 
i and i^) and belongs to ^"^((0, T) x fi) for all q G [1, 00). Thus, for the sake of applying Moser iterations, 
relation (j4.28p shows that we can proceed by working exactly as in the linear case (with forcing term 
given here by ^iJ^/^). Actually, by (HT5)) . (IT^U)) and interpolation we have cj)§^^'^ G L''((0,r) x Q) for a 
suitable r > 3. Thus, applying the results in |17j . we obtain the bound 

l|i^l|L~((o,T)xa) < c{N, m, e, 6), (4.29) 

and this estimate depends on the approximation parameters, but is uniform with respect to time. Finally, 
applying the theory developed, e.g., in |16j Section 4], we obtain 

ll'^llc°([o,T]xn) < c{a,N, 111,6,6), for some a > 0, (4.30) 

globally in time. 



4.4.5 Estimates based on energy dissipation 

The energy bounds analogous to (|3.7p . (I3.10p can be deduced from (|4.16p exactly as in Subsection [ 
Actually, although Q may not lie in L°° at this stage, the term 



C 



§(V,u, 



does not suffer any summability loss since we can suppose that G has zero gradient when 
Hence, we get 

|VxU||L2((o,T)xn;_R3x3) < C, (5^^''||Va:U||ir((o,T)xn;_R3x3) < C, 



<c. 



L2((o,T)xO;_R^;„%) 



where all constants are independent of the parameters m, N, e, and 6. 
By (|4.1ip . (|4.33p and interpolation, we also infer 

ll^llLiO((O.T)xO;i?,fy><^o)) ^ ll^a:Q|!Lio/3((o,T)xn:i?27)) 

Analogously, from (I4.10p and the first (|4.3ip we get 

l|u|!L"'/3((0,T)xn;_R3)) < C. 



< c. 



is large. 
(4.31) 
(4.32) 

(4.33) 

(4.34) 
(4.35) 
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Coupling (|4.34p . (|4.35p . and (|4.3ip . we obtain a uniform bound for the coupling terms in (|4.4p . namely 

||u. V.Q-§(V.u,Q)||^,„,,((^_^^^^^^3><^3^^^ < c, (4.36) 

whence (|4.32f[4?33|) and a comparison of terms in (|4.4[) lead to 

\\dtQ\\Lio/7^^o,T)xn;Rl-l,) < c. (4.37) 
We notice however that estimates (|4.34H4.37|) will be improved in the sequel. 

5 Convergence to the limit system 

The bounds derived in the previous section, being global in time, are sufficient for extending the approx- 
imate solutions to the desired existence interval. For this reason, and in order to avoid technicalities, we 
shall directly assume that solutions are defined over the whole (0,T) already at the approximate level. 
Our ultimate goal is to perform the limits 

m — >■ (X), N — )■ oo, e — >■ 0, and, finally, J — > 0. 
5.1 Uniform bounds on the temperature 

Testing (|4.16|) by —(1 — i?)^", a > 0, it is not difficult to arrive at the analogue of p.lOp . namely 

VM + < c, for any a > 0, (5.1) 

Indeed, all terms on the right hand side of (I4.16P are uniformly bounded in L^. In particular, C 

S(Va;U, Q) has the same regularity of V^u because we can suppose G to be have zero gradient in the set 
where |Q| is large. Note also that regularity and strict positivity of the (approximate) d are essential in 
order that this procedure makes sense. 

Coupling the information coming from (I4.12p and (|5.ip . which are both independent of all approxi- 
mation parameters, and using interpolation arguments together with the Sobolev imbedding W^'^(i7) ^ 
i^(n), we arrive at 

II^I|li((o,t)xO) < c for any 1 < g < 5/3. (5.2) 
Moreover, another application of interpolation (see also [SJ Sec. 4]) permits to obtain 

l|Va;t?|!L9((o,T)xn;fl3) < c for any 1 < g < 5/4. (5.3) 

Next, thanks to the lower bound established in (|4.24[) and to the Holder regularity (|4.30|) . we can (rigor- 
ously) divide the heat equation (|4.7p by i) to obtain the entropy relation 



d, 



( log(i^) + U's{d)G{Q)) + u • V, ( log(i9) + U's{d)G{Q)) - div, (^^-^) ^^'^^ 
= I I V,u + V^u|' + ^|V,ur + [r(^)],|H„,,|2 + ^iMjv.^p^ , 
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which is still an equality at this level. 
Integrating (|5.4[) . we deduce that 



ess sup II log(i?)||Li(n) < c, II log(i?)I|i2(o,T;M^i.2(n)) < c. 

tG(0,T) 



(5.5) 



Finally, we need an estimate on the time derivative of d. To get this, we go back to (|4.7p (or, equivalently, 
multiply (|5.4[) by d) which we rewrite in the form 



9t(i? + M5(i?)G(Q)) +u-\7,(d + Msi^)G{Q)^ - div, (k(i9)V,i9) 
= (AW) - U^{^)^) {dtG{Q) + u • V,G( 



(5.6) 



V,u + V^ul +<5|v,ur + [r(i?)] 



where we have set 



Us{r)r dr. 



Then, testing (|5.6p by a test function (f) S If''(ri) and using the previous estimates (|4.31H4.35|) and (|5.2p . 
and the fact that, by (fLBHTT)) . Ms{'d) ~ t?^/^ for large i9, it is easy to check that 



<c, 



(5.7) 



uniformly with respect to all approximation parameters. Since G and are smooth and nonnegative 
and 9tQ is already estimated in (|4.37p . we finally obtain 



\\dt'd\\L^{a,T:H-^(n)) < c- 



(5.8) 



5.2 The limit m — t- oo 



It is convenient to start with the limit for m -t- oo, with the other parameters iV, £, and 5 fixed. In such 
a way, the velocity field remains regular at this stage, with facilitates the limit passage considerably. As 
a result of this step, we obtain that /(Q) lies in L°°(0, T; L^(ri)) in the limit. This yields, in particular, 
uniform boundedness of Q by a constant independent of the parameters N , e, and 5. 
We start by recalling that, in accordance with (|4.9p . 



(5.9) 



ess sup / /„i(Qm) da; < c, 

*G(0,T) Jn 



whereas, by virtue of the standard identity (see for instance [24] ) 



C 



dfrr 



fn 



dfr. 



combined with ((432|) and (|434l) . we infer that 

l./m(Qm)llLV3((o,T)xa) < C. 



(5.10) 
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We are ready to perform the limit for m — )■ oo. To this end, we denote by {u^, i?m}m>o the family 
of approximate solutions. First, by virtue of (|5.9p and (|5.10[) . we have 

/r«(Q„J ^7(0) weakly-(*) in {0 , T ; L^Q)) and weakly in L5/3((o,r) x f]), 

at least for suitable subsequences. 

Next, using (|4lT1) . (|03l) . (g^H), (|437ll . and the Aubin-Lions lemma (cf. [ll[18]), we get 

Qmit, •) ^ Q{t, .) in, say, L^{n; R^^^^ ,,) for a.e. t e [0, T], 
VxQ™ ^ V,Q (strongly) in L«((0,T) x n;R^^) for aU q e [1,10/3), 

and 

AQ„ ^ AQ weakly in ^^((o, T) x R^^^l^). 
Our ultimate goal is to show that 



,x) e V[f] for a.a. {t,x) G (0,T) x 17, (5.11) 
yielding, in particular, the desired conclusion 



Using (|5.10p . the monotonicity of the sequence fm, the pointwise convergence of {Qm}m>o, and 
Fatou's lemma, we obtain 

ess sup / /mo('Q) dx < c for any fixed mo, with c independent of toq; 

tG(0,T) 

whence, by the Levi theorem applied for toq ~^ oo, we conclude that 



ess sup / f{Q) dx < c, 
te{o,T) Jn 

yielding ([5lT|l . 

Let us notice that, as a consequence of (|5.1ip . we have in particular 

Recalling (|4.33p and applying the standard Gagliardo-Nirenberg inequality (cf. fSP, p. 125]), we then also 
have 

llV:rQ||i4((o^T)xf2;fl27) < C, (5.13) 

uniformly in N, 5 and e. As a consequence, we can also improve (|4.36l - l4737|) as follows: 

||u.V.Q-S(V,u,Q)||^,„/,,^(^^^)^^^^3><^^^ <c, (5.14) 
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5.2.1 Strong convergence of the temperature and the Umit system 

Using (|5.2H5.3p . (|5.8p . and the Aubin-Lions lemma (cf. dUIH]), we deduce that 

z9™ ^ ?? in i«((0,r) X n) for a certain q > 1. 

Then, the uniform estimates derived above permit to take the hmit rn — )■ oo in the momentum equation 
to obtain 

d f 

— / u • V dx (5.16) 
Jn 

= / [u]s (g) u : Va;V dx-6 : Va;V dx ~ ^(i?)(V;j;U + V^u) : Va;V da; 



2^ (Ms : Q) 
for any v G Xn, where 



'0 V, 



VrV da; 



Va;V dx, 



dG{Q) 



Analogously, we can take the limit m — ^ oo in the director equation, obtaining 

dtQ + (u • V,)Q - S(V,u, Q) - [r(?9)]eH5. 



(5.17) 



(5.18) 



Actually, we remark that u is still a smooth vector field in the limit. On the other hand, we cannot prove 
strong convergence for H and Vx"^- For this reason, the heat equation cannot pass to the limit m ~^ oo 
in the form (j4.7p . For this reason, it is convenient to replace it by a partial form of the energy balance 
at this stage. 

To do this, for fixed m we multiply (14. 4[) by ~Mm,5 and sum the result to (|4.7p . Note that we 
do not sum the energy contribution coming fi'om the momentum equations at this level. A number of 
integrations by parts similar to those performed in Section 13.31 then permit to deduce 



dt ( i|V 



(5.19) 



-diV:r ( U„, ( -|Vj;(I 



iv.(K(i9™)V,An) -div,(v,Q,„ :S(V,u„,(Q„j) - div,([r(i?™)],V,(l 



-§(V,U„,Q,„) 



( 







: V^u„ 



Vj;U™ + V^u„ 



Now it is possible to take the limit m — ^ oo in the above relation. Actually the worst terms are the 
quadratic ones on the right hand side. However, they do pass to the limit since the velocity still takes 
values in the finite dimensional space ATjv at this level (and, hence, it is a smooth function). We then get 



dt 



1 



Us{i))-m's{i)))G{Q) + ^ 



(5.20) 
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+diVa; 



-div,(K(i9)V,i9) -div,(v,Q:S(V,u,Q)) - div, ([r(i?)],V,Q : H,-) 



-§(V,u, Q) 



Due to the lack of strong convergence for H and Vj^i?, also the m-limit of the entropy equation (|5.4p has 
to be written in the form of an inequality. Indeed, a standard semicontinuity argument yields 



9t(log(i?) + U'smCiQ)) + u • V,(log(i9) + U'sid)G{Q) 
1 /AiW 



— diVa; 



(5.21) 



> 



v,u + v>r + (5|V,ur + [r(i9)]e|H|- 



5.3 The limit N ^ oo 



Our next goal is to let the Galerkin parameter N — oo. In this step, the most difficult points consist 
in taking the limit of (|5.16p . in order to let it converge to the Navier-Stokes system, and of the internal 
energy equation (j5.20p . To achieve these limits, we need to prove a strong convergence for V^un- With 
this purpose, we denote {ujv, Q^r, i?Ar}Ar>o the family of approximate solutions and recall that 



^JvllLoo((o_T)xO;_R-'y><,^ (,) - ^' 



(5.22) 



with c independent of N, e, and 6. Then, we take v = ujv in (|5.16p and correspondingly test ()5.18|) by 
Hat. Repeating the usual cancellations, we then arrive at 



i|ujvp + i|V,QArp + /(Qjv) ) dx~ I Us{dN)C 



N 



IN 



(atQftr + uwV.Qjv) da; (5.23) 



5\V^\in\ 



/i(^?iv) 



V,UAr + V^uw|' + [r(i?w)],|Hwn d.T = 



where we have used the identities (|1.33p and 

(UAT • V^Qat) : AQat + (V^Qjv V^Qjv) : V^uat dx = 0. 



Jn 



Next, we integrate (|5.23p in time and compute its supremum limit as — ?► oo. We obtain 



lim sup 

jv^-cx) Jo Jn 



lim 

N-yoo 



^ ' ' r , M(^iv) 



(5|Va;UAr| 



Jn 



V,UAr + V>Ar| +[T{dN)]emN\' j dx dt 
• (^tQAT + U^v • V,Q^) dx dt 



(5.24) 



(5.25) 



UN 



■ lim inf 

N^oo 



i|uAr(T)p + i|V,QA.(T)p + /(Q^) ) dx + 



lli^oW + l\VMs\' + fiQo)) dx 
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n 



{dtQ + u • V^Q) dx dt 



Actually, the existence of the limit in the second line above is ensured by the pointwise convergence of 
•On, hy the boundedness of Us and G together with their first derivatives, and by estimates (|5.2p and 
(|5.14l45TT5)) . Moreover, the liminf on the third line above can be estimated by using semicontinuity of 
norms with respect to weak convergence. 

Now, we observe that it is possible to take the limit iV — ^ oo in (|5.18|) . Analogously, we can write the 
limit of (|5.16p : however, the second integral on the right hand side has to be temporarily written as 



iVjjul'-^v^u : Vj;V dx (5.26) 

in the limit, since strong convergence of V^^u is not achieved yet. 

Then, we take v = u in the iV-limit of ((5?T6l) and test the TV-limit of (fSTTS]) by H. Note that this 
procedure is rigorous (and gives rise to an equality) thanks to the fact that we have that V^^u is bounded 
in L'', where r > 3, uniformly in N and, consequently, the coupling term 

u- V,Q-§(V,u,Q) 

in (|5.18p lies in even in the limit. Hence, equation (|5.18[) can still be read as a relation in and use 
of the L^-test function HI is consequently permitted (note that this will no longer be true in the limit 
5^0). This is exactly the reason why the r-Laplacean regularization has been added in the momentum 
equation. 

This procedure permits to achieve, in the limit — !■ cxd, the analogue of (j5.23p . Comparing with 
(|5.25|) we then obtain 



(5.27) 



limsup r / f<5|V,U^r + ^^^|V,UAr + V>Ar|' + [r(z?Ar)],|HAr|2) dx dt 

<^ J (^(JlV.ul'-V.u . V.u + ^ I V,u + V^u|' + [r(i9)],|H|2^ dx dt, 

whence a simple monotonicity argument permits to deduce the strong convergence 

V^ua, ^ V^u in L''((0,T) X r2;i?9). (5.28) 

Using this relation together with (|5.13p . (|5.15p . and the Aubin-Lions lemma (cf. |T1[TH]), we can also take 
the limit A — cx) of (|5.20p (and in particular of the quadratic terms on the right hand side). 

Moreover, (|5.28p permits to identify the limit of | V:cu|''~^Va;U in the momentum equation (cf. (I5.26p ). 
In particular, the Navier-Stokes system (|5.16p passes to the desired limit as A cx). 

Next, we notice that the (infimum) limit N oo can be taken in the entropy inequality ()5.2ip as in 
the previous section. Thus, to complete the passage to the limit w.r.t. A^ — oo it is sufficient to recover 
the total energy balance (cf. (|1.37p ). With this aim, we can notice that, in the limit A -> oo, (|5.16p is 
no longer a relation in the finite-dimensional space Xjy, but can rather be interpreted as a true PDE, (at 
least) in the distributional sense. In other words, we can rewrite it in the form 

dtu + div.j,{[u]s u) = div^cr, (5.29) 
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where 



(7 = At(i?) (V^u + V^u) +5|V^ur "^V^u-pl 



(5.30) 



Then, multiplying (15. 29^ by u and adding the result to (|5.20p . standard integrations by parts permit to 
get the total energy balance ()1.37|) (where of course, at this level, the stress a still contains the regularizing 
contribution (5| V2:u|'"^^Vxu). 



Remark 5.1. In order to achieve the total energy balance (|1.40p it is crucial to remove first the Galerkin 
approximation ( otherwise the kinetic energy contribution is projected on the finite- dimensional space Xn ). 
For this reason we need that, after taking the limit N oo, we still have sufficient regularity to use u 
as a test function in the momentum equation in order to get a kinetic energy equality. This regularity is 
properly provided by the additional term 5| Va;u|''~^Va;U. 



5.4 The limits e and 5^0 

After taking the limit — )> cxd, we still have to remove the regularizations coming from the parameters 
e (appearing in the coefhcient [r('(9)]e in (|4.4p ). and S (appearing in the mollification of the initial data, 
in the function Us in the approximation of H, and in the regularizing terms [u]^ ® u and (5| Va;u|''~^Va:U 
in the momentum equation). 

Since the total energy balance (|1.40[) is already achieved and no quadratic terms are present in its 
right hand side, letting e ^ and subsequently 5^0 does not give rise to any additional difficulty. 
Actually, most of the argument can be carried out just by adapting the procedure, based on the uniform 
a-priori bounds of Scction[3l used before to let m, N ^ oo. The only point which requires some additional 
care is letting (5 — in the momentum equation. However, this procedure is absolutely analogous to the 
argument outlined in [8j Sec. 5.2] to which we refer the reader for details. We just note that we need 
here the restriction r < 10/3 (cf. Section |4?T|) on the exponent of the additional viscosity term (up to this 
point we only used that r > 3). 
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